CMB bounds on dark matter annihilation: nucleon energy-losses after recombination 
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We consider the propagation and energy losses of protons and anti-protons produced by dark 
matter annihilation at redshifts 100 < z 2000. In the case of dark matter annihilations into 
quarks, gluons and weak gauge bosons, protons and anti-protons carry about 20% of the energy 
injected into e* and 7's, but their interactions are normally neglected when deriving CMB bounds 
from altered recombination histories. Here, we follow numerically the energy-loss history of typi- 
cal protons/antiprotons in the cosmological medium. We show that about half of their energy is 
channeled into photons and e*, and we present a simple prescription to estimate the corresponding 
strengthening of the CMB bounds on the dark matter annihilation cross section. 



PACS numbers: 95.30.Cq, 95.35.-|-d, 98.80.Es 



I. INTRODUCTION 

Astrophysical and cosmological observations provide 
compelling evidence that about 85% of all the matter 
in the Universe is in the form of Dark Matter (DM) , an 
elusive substance which is currently searched for with a 
variety of observational and experimental channels at col- 
liders, in underground detectors, or via indirect signals 
from DM annihilation or decay [l[ . 

Cosmic microwave background (CMB) anisotropy and 
polarization data provide interesting constraints on the 
properties of DM particles Secondary particles 

injected via DM annihilation (or decay) after recom- 
bination, around redshift z ~ 0(600), would in fact 
inevitably affect the recombination history of the Uni- 
verse and widen the surface of last scattering, which is 
severely constrained by CMB observations as discussed 
in Refs. [^-lll: and more recently in Refs. [9l-[T5|. 

One of the main reasons of interest for these constraints 
is that, in contrast with other indirect searches, they do 
not rely on knowledge of astrophysical DM structures, 
affected by the complex aspects of non-linear gravity as 
well as complicated feedback due to baryons. The CMB 
probe is thus as reliable as the description of the basic 
atomic and nuclear/particle physics processes involved 
is. Its robustness (and astrophysical independence) mo- 
tivates further efforts to assess and improve the error 
budget of bounds obtainable with this technique. 

It is typically assumed that protons are highly pene- 
trating and poor at transferring energy to the intergalac- 
tic medium (IGM) — a misnomer, since no Galaxies have 
formed at such high redshifts — and their energy release 
to the medium is neglected 0, [Ij. In this article, we 
estimate the additional energy released to the gas by the 
interactions of the high-energy protons and antiprotons 
formed among DM annihilation final states. 

This article is structured as follows: In Sec. [Ill we 
review the basic physics substantiating the two points 
above. In Sec. IIIII we describe our computational tech- 
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nique and present our results in Sec.|IVl Finally, in Sec.lVl 
we conclude. 



II. PHYSICAL PROCESSES 

Protons and antiprotons carry a significant fraction 
of the overall energy emitted in DM annihilations into 
quarks, gluons or weak gauge bosons. Typically, this 
amounts to ~ 20% of the energy channelled into e* and 
7's, see for example Fig. 4 in ^M,^ and a fraction of this 
energy will be inevitably transferred to the IGM. Neu- 
trons and antineutrons decay very fast and behave practi- 
cally like protons and antiprotons, while (anti-)deutcrons 
and heavier nuclei are produced in negligible amounts in 
DM annihilations. 

At the epochs of interest here, which correspond to a 
redshift z = 0(10'^), p/p's propagate in a medium which 
is eight or nine orders of magnitude more dense that at 
present epoch, with typical proton and Helium gas densi- 
ties up to 0(10^) particles/cm"^. Even neglecting the in- 
teraction with photon baths of densities of C'(10^^)cm^'^, 
a typical p-p inelastic cross-section of 30 mb yields a col- 
lision timcscale lower than the age of the Universe at 
decoupling. Hence, we expect that in general the prob- 
ability for a nucleon to interact within a Hubble time 
is large, and that a significant energy deposition takes 
place. 

More specifically, p/j5's will undergo the following pro- 
cesses: 

• Coulomb scattering of p/p on IGM electrons [TtI "^ 

• Thomson scattering of p/p on IGM photons [l^ 

• Inelastic scattering of p/p on IGM protons [loj 



^ As lower cutoff in the Coulomb logarithm wc adopt 10 eV. 
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FIG. 1: Upper panel: Total inelastic hadronic cross-sections 
in lab frame as function of the p and p kinetic energies. 
Where applicable, we also show the annihilation part of the 
total cross-section separately. Lower panel: Fractional en- 
ergy losses relative to Hubble rate, at redshift z = 1000. The 
non-annihilating and annihilating parts of the hadronic cross- 
sections are shown separately. 



• Inelastic scattering of p/p on IGM Helium [T^. 

In Fig. [1] we show the inelastic scattering cross-sections 
as a function of the p/p kinetic energy (top panel), and 
the corresponding fractional energy loss rate E~^dE/dt 
times the Hubble time at redshift z = 1000 as 

function of the p/p kinetic energy (bottom panel). Val- 
ues of order unity or even larger indicate potentially 
large effects. In case of anti-protons, we plot annihi- 
lating and non-annihilating rates separately. For non- 
annihilating inelastic interactions, the final energy dis- 
tribution of p/p is taken to be constant in the physi- 
cally accessible regime rup. . . Ep {Ep is the energy of the 
primary particle, rUp the p/p mass). Consequently, the 
mean energy loss during a non-annihilating collision is 
(AE) ~ O.STp, with Tp = Ep — rUp being the initial ki- 
netic energy. The energy released during an annihilation 
event is AE = E„ + ■ 



losses are dominated by inelastic p-p and p-p scattering 
without annihilation; at lower energies p-p annihilation 
and Coulomb losses become relevant and dominate for 
non-relativistic particles. 

In the next section we shall follow the evolution of a 
nucleon pair from a hypothetical toy annihilation process 
XX ~^ PPi computing the energy loss of the daughter nu- 
cleons down to relatively low-rcdshifts of a few hundreds. 
Obviously, only part of this energy will be absorbed by 
the gas. Following in detail the energy degradation down 
to the energy-transfer processes to the gas — heating 
and ionization — goes beyond our purposes (it is worth 
mentioning that several efforts are being put into a more 
realistic treatment of the related physics for the e^ and 
7's, see e.g. H^])- Rather, we will content ourselves with 
estimating the energy that is released into energetic elec- 
trons and photons (essentially, the fraction of the energy 
lost to stable particles other than neutrinos). We will 
express this as the electromagnetic fraction fe±-y of the 
energy initially injected as nucleons, which we will pre- 
cisely define below. This approach is exact as long as the 
CMB bounds are "calorimctric" and the precise form of 
the spectrum of electromagnetic particles is irrelevant in 
determining the ultimate fate of the energy injected. 

Previous investigations suggest that this should pro- 
vide a reasonable first approximation to the true result; 
for our purposes the accuracy is expected to be at the 
level of - 30% (see Fig. 2 of [H, z - 600). Since we 
are already concerned with a correction to the basic re- 
sults, we shall adopt this Ansatz which greatly simplifies 
the problem. The only potentially problematic case is the 
one of Coulomb reactions of the (anti-)protons, which can 
produce low-energy electrons (which may behave differ- 
ently from high-energy ones). However, the legitimacy 
of our approximation is supported by two arguments: i) 
this energy-loss rate has only a logarithmic dependence 
on the minimum kinetic energy; ii) it only matters for 
non-rclativistic (anti-) protons, for which most of the en- 
ergy coming from electrowcak-scale WIMPs has already 
been dissipated. In fact, with the exception of very light 
(GeV scale) DM particles, protons are typically born rel- 
ativistic or semi-relativistic. 



III. COMPUTATION 

The energy loss processes sketched above are of two 
distinct classes: continuous energy losses (with small en- 
ergy transfer per collision) when protons interact with 
electrons or photons; or catastrophic losses, when they 
undergo an inelastic collision or an annihilation (for an- 
tiprotons). For the moment we shall ignore elastic colli- 
sions, we will comment on their role at the end of Sec. lIVI 

Continuous energy losses in the range of our interest, 
100 < z ^ 2000, can be described via the following dif- 
ferential equation 



At energies > 10 GeV, the 
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The last term at the RHS describes adiabatic energy- 
losses in the expanding universe, the other two terms 
describe Coulomb and Thomson energy-losses. In the 
present work, Thomson losses will be neglected, which 
from Fig. [T] is clearly justified at the energies of interest. 

Hadronic processes arc included on top of the contin- 
uous energy losses with a Monte Carlo (MC) simulation. 
The MC starts with a fixed p/p energy Einj at some ini- 
tial redshift zinit and tracks the history of the particle as 
it moves to lower redshifts. The redshift of an individual 
scattering process is inferred by sampling from a survival 
function, which is determined by solving a differential 
equation whose derivate is given by the inelastic scatter- 
ing rates. The amount of electromagnetic and hadronic 
energy lost by each particle as function of redshift z is 
recorded. 

At any given redshift z, a fraction fe±i{z) of the 
hadronically injected energy is actually channeled into 
electromagnetic form (energetic and photons, see dis- 
cussion above). More specifically, for annihilating DM, 
/e±^ is defined as 

dE 

ee±7(^) = ^ = 2A/;,/,±^(z) • . (1 + zf , (2) 

where e^±^{z) denotes here the energy released into elec- 
tromagnetic form per comoving volume and per unit 
time, and n~^fl is today's number density of DM parti- 
cles. This fe±-yiz) can be derived from the integral 

feijiz) = dz' g,±^{z, z') j^^^,)(-^^^-)3 ' 

where g^±y{z,z') denotes the fraction of the energy of 
a single DM annihilation event at redshift z' that is re- 
leased into electromagnetic energy at redshift z per unit 
time. We will take Zmax = 10000 throughout, such that 
ee±^(z) at z < 1000 is independent of ^max- 

We split ff^±^{z), and analogously the function ge±^(z), 
into three parts according to 

/e±7 = /C ou. I Q^had (/inol. +/a„n.) , (4) 

with the three addenda at the RHS corresponding to 
Coulomb losses, inelastic scattering and annihilation, re- 
spectively. These /'s denote the total (anti-)proton en- 
ergy lost to secondary particles via the scatterings. All 
of the energy lost by Coulomb scattering, and a fraction 
Q^had of the energy lost during inelastic scattering and an- 
nihilation is channeled into electromagnetic energy and 
hence contributes to fe±f 

A reasonable value is ahad — 0.5. After a hadronic 
scattering, the energy eventually carried by and 7 
contributes fully to fe±-y, the energy carried by i^'s is 
completely lost, and only a fraction of the energy carried 
by nucleons will be transfered into 6^7 at a later stage. 
We will neglect this latter contribution to /e±-y, since 
it is further suppressed by the small fraction of energy 
released into nucleons during inelastic scattering. 



As a proxy for the energy released in the different chan- 
nels we can take the gluon-gluon (or charmed quark) 
channels of Fig. 4 in [1^, where about 46% of the en- 
ergy is found to be channelled into and 7, with 
^ 27% of this amount going into nucleons. We also 
checked that these figures do not change appreciably as 
a function of M-^ (M. Cirelli, private communication). 
This is also consistent with the expectations from non- 
relativistic proton-antiproton annihilation, see e.g. Table 
2 in [2l| . A very simple justification of this value can be 
obtained in the limit where annihilation final states are 
just made of pious, with isospin-blind production yields 
(equal quantities of 7r°, 7r+, tt"). One would expect full 
channeling of 7r° energy into e^j and about 1 /6 of the en- 
ergy stored in , hence a weighted average of 44%. We 
thus believe that ahad — 0.5 is a reasonable benchmark, 
probably affected by a 20% uncertainty. While this may 
be a crude description of the actual process, it is enough 
for providing a first estimate of the hadronic correction 
to CMB bounds on annihilating dark matter. 

On a more technical note, to obtain gcou.(^,^')j 
5inei.(-Zj z') and ^annX-^i z'), wB simulate particle injection 
at redshifts z' = 30-10000 and record what fraction of the 
energy is released in certain redshift intervals zq . . . zi^ di- 
vided by the corresponding time interval At = to ~ ti. 
From these g{z, z'Ys, we can then derive the ff,±^{z) from 
Eqs. Q and gl). 



IV. RESULTS AND DISCUSSION 

In Fig. [2] we show the main results of this paper. In the 
top panel, the resulting different components of the frac- 
tional energy loss of (anti-)protons are shown as function 
of their injection energy at redshift z ^ 600. Notably, 
the annihilation of anti-protons releases also the energy 
of the target proton, which yields values /ann. > 1 when 
considering the anti-proton channel alone. The bottom 
panel shows the fraction of the injection energy that is 
channeld into electromagnetic energy, /e±^, as function of 
redshift z and for different kinetic injection energies Tp/p. 
For energies Tp/p < 10 GeV (Tp/p > 10 GeV) about 
fei-y ^ 50% (40%) of the energy fraction channelled into 
nucleons will eventually by released as electromagnetic 
energy at the most relevant redshifts z = 500-700. 

Most importantly, for DM masses below ^ 100 GeV, 
the /e±^(z)'s are approximately independent of redshift. 
This implies that the released electromagnetic energy is 
simply a constant correction to the prompt electromag- 
netic energy directly released during DM annihilation. 
In the current literature (see e.g. Ref. Q), the fraction 
of the total injected energy finally deposited in the gas 
as heat, ionization and excitation energy is denoted by 
/(z), and determined neglecting annihilation channels 
into hadronic final states. As discussed above, CMB 
bounds are to a good approximation "calorimetric" , and 
the precise form of the spectrum of electromagnetic par- 
ticles can be neglected at first order in determining the 
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FIG. 2: Upper panel: Fractional energy losses at redshift 
z ~ 600, for different kinetic injection energies Tp/p. We show 
the contributions from annihilating and non-annihilating 
hadronic scatterings and from Coulomb losses separately. 
Lower panel: Electromagnetic (i.e. energetic and photons) 
energy released during p/p propagation, as function of red- 
shift z, for the toy process xx ~^ PP with different dark matter 
masses — 1.0 GeV to — 1 TeV. In all cases we show 
Monte Carlo results (crosses) as well as interpolations with 
splines of degree two (lines) . 



ultimate fate of the energy injected. We can then give a 
simple recipe to correct these /'s: 



/(2)->/(z)(l + /e±^- 



7+e 



(5) 



where i?NN ^■^d E^_^_f.± denote the prompt energy in nu- 
cleons and in electromagnetic energy, respectively, pro- 
duced during an annihilation. As shown in the bottom 
panel of Fig. [2l we find that /e±7 — 0.5 for injected 
p/p energies well below ~ 10 GeV, i.e. for most DM 
masses of interest. 

Throughout, we neglected elastic p/p scattering, which 
is well justified for our purposes: At energies below Tp < 



1 GeV, elastic p-p scattering will only redistribute energy 
that finally is lost via Coulomb scattering. At energies 
above Tp > 1 GeV, elastic p-p scattering is subdominant 
50% of the total p-p cross-section); even rare collisions 
with maximal energy transfer do not alter the deposition 
history drastically since the energy dependence shown in 
Fig. [2] is weak. In case of p-p scattering, the elastic cross 
section is smaller than 40% at all energies; the main effect 
of collisions with large energy transfer would be a slight 
increase in anti-proton annihilation rate. 



V. CONCLUSIONS 

It has been recently recognized that non-standard 
sources of heating and ionization of the IGM can be 
probed via CMB anisotropics. This effect has been used 
to derive bounds on DM annihilation cross-section and 
decay lifetime. To a very good approximation, these 
bounds do not depend on "astrophysics" (structure for- 
mation, DM halo structure, star formation etc.). The 
nature of this probe, whose robustness relies essentially 
on the accuracy of the description of atomic and nu- 
clear/particle physics, motivates efforts to assess and im- 
prove the error budget of bounds obtainable with this 
technique. 

In astroparticle physics, energy losses of high energy 
nucleons propagating in the intergalactic medium are 
usually neglected outside the realm of ultra-high energy 
cosmic rays. Energetic nucleons propagating in a high- 
z cosmological medium over Hubble times, on the other 
hand, do suffer significant energy losses. A physically in- 
teresting example is provided by nucleon by-products of 
dark matter annihilation, which can release a significant 
fraction of their energy at high-z (~ 0(10'^)). 

The aim of this paper was to provide a first estimate 
of this process, which goes in the direction of improving 
the CMB bounds for a number of channels. Under rea- 
sonable approximations, we found that the Eq. ([5|) is a 
fair description of the main effect, for a value /e±-y — 0.5 
over a quite large parameter space. We expect thus that 
for DM annihilation channels into gluons, gauge bosons 



and quarks for which 



0.2, this should change 



the CMB bounds on DM at the 10% level and should be 
included in forthcoming Planck data analyses. 
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